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Abstract. 

A study of the luminosity function of 36 Abell clusters of 
galaxies has been carried out using photographic plates ob- 
tained with the Palomar 1.2 m Schmidt telescope. The relation 
between the magnitude Mi of the brightest cluster member 
and the Schechter function parameter M* has been analyzed. 
A positive correlation between M* and Mi is found. However 
clusters appear segregated in the Mi-M* plane according to 
their Rood fc Sastry class in such a way that on average Mi 
becomes brighter while M* becomes fainter going from late to 
early Rood fc Sastry and also Bautz fc Morgan classes. Also 
a partial correlation analysis involving the magnitude Mio of 
the 10th brightest galaxy, shows a negative intrinsic correlation 
between Mi and M* . These results agree with the cannibal- 
ism model for the formation of brightest cluster members, and 
provide new constraints for theories of cluster formation and 
evolution. 

Key words: Galaxies: clusters of - galaxies: evolution - 
galaxies: photometry - galaxies: luminosity function - galaxies: 
statistics 



1. Introduction 

The observed luminosity function (LF) of galaxies provides 
the most important means to test theories of galaxy formation 
and evolution. According to a reformulation and extension of 
the Press fc Schechter (1974) theory presented by Bond et al. 
(1991), the direct hierarchical clustering of primordial overden- 
sities is responsible for the onset and evolution of a self-similar 
mass distribution leading to the observed Schechter-like galaxy 
luminosity function. 

Moulding of the luminosity function by galaxy merging 
has recently been exploited (Cavaliere fc Menci 1993 and refs. 
therein) to reconcile the local LF (Efstathiou, Ellis fc Pe- 
terson 1988; Binggeli, Sandage fc Tamman 1988) with faint 
galaxy counts (Tyson fc Seitzer 1988) and redshift distribu- 
tion (Broadhurst, Ellis fc Shanks 1988). 

Galaxy merging, or galactic cannibalism, was originally in- 
troduced by Ostriker and Tremaine (1975) as a mechanism 
of growth of brightest cluster members at the expense of the 



other massive galaxies, which are most affected by dynami- 
cal friction. According to the model of Hausman and Ostriker 
(1978), as the evolution of the cluster proceeds this selective 
depletion should push to lower luminosities the turnover point 
between the steep high-luminosity falloff and the flatter faint 
end of the LF. This should cause a negative correlation be- 
tween the magnitude Mi of the brightest cluster member and 
the characteristic magnitude M* of a fitting Schechter (1976) 
LF. 

From a study of 12 rich clusters, Dressier (1978) derived an 
indication that Mi and M* are anticorrelated. However, in a 
subsequent study of 9 Abell clusters, Lugger (1986) (L86) con- 
cluded that the correlation is not statistically significant and, 
in a more recent study of 12 Abell clusters, Oegerle and Hoessel 
(1989) (OH89) find no evidence for any relation between Mi 
and M* . These results are taken, by the respective authors, 
as indications against the effect found by Dressier (1978). It 
should be noted that in both these studies a positive Mi-M* 
correlation is found, though it is not statistically significant 
due to the small number of clusters analyzed. 

Since the properties of individual clusters are strongly af- 
fected by random fluctuations, they become meaningful only 
when defined statistically using a large number of clusters. 
Thus a project has been undertaken for studying in a uni- 
form manner a large sample of nearby galaxy clusters (Flin 
et al. 1988, Trevese et al. 1992 (T92), Flin et al. 1995), and 
deriving properties such as number density profiles, morphol- 
ogy, galaxy orientations, luminosity functions (LF) and their 
possible statistical relations. 

The results concerning galaxy orientations in a sample of 55 
clusters (Trevese, Cirimele and Flin 1992) indicate the special 
role played by first ranked galaxies during the evolution of 
clusters. 

In the present paper we restrict our attention to the spe- 
cific problem of the Mi-M* relation and present the results 
obtained from a sample of 36 Abell clusters, more than three 
times larger than each of the previous samples. From this we 
obtain a statistically significant evidence of a new type of neg- 
ative Mi-M* partial correlation, related to the fact that Mi 
becomes brighter and M* fainter in going from irregular to cD 
Rood fc Sastry (1971) (RS) cluster types or from type III to 
type I in the Bautz fc Morgan (1970) (BM) classification. 
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Our data are derived from 10-inch photographic plates 
taken by P. Hickson with the 48-inch Palomar Schmidt Tele- 
scope to analyze a sample of 64 Abell clusters. The fields were 
selected according to the criteria specified in Hickson (1977) to- 
gether with details about the emulsions and filters employed. 
The resulting photometry corresponds to the red F-band of 
Oemler (1974). The sample was not statistically complete but 
it was designed to cover all cluster morphological types, to 
perform statistical studies of each cluster type. Other clusters 
appearing in the same plates, some of which are not Abell 
clusters but belong to the Zwicky catalog, were also added to 
the sample, thereby reaching a total of more than 100 clus- 
ters in all. Plates were scanned with a PDS 1010G in Rome, 
with pixel sizes ranging from 10 to 25 (im depending on the 
cluster distance. Automatic identification of objects and their 
classification as point-like or diffuse are described in T92. Total 
magnitudes are computed from the flux integrated in a circular 
aperture whose radius is Ri = 1.5ri, where n is the first mo- 
ment of the intensity distribution (see T92). The magnitude 
defined in this way corresponds on average to an isophotal 
magnitude at 24 mag ■ arcsec~ 2 , with the advantage that n is 
less noisy than the corresponding isophotal radius. The signal 
to noise ratio is S/N ^ 100 for a few objects brighter than 
F = 12 mag, about 25 for F PS 14 mag and falls to about 5 
for F ^ 18 mag. 

Relative photometry has been obtained for 55 of the above 
clusters (Trevese, Cirimele and Flin 1992), while the zero of 
the magnitude scale has been established for 36 clusters, using 



published photometric data. For 27 out of 36 clusters we used 
r band data from Hoessel, Gunn fc Thuan (1980) and Hoessel 
fc Schneider (1985); for 8 of the remaining clusters we used V 
data while for A2593 we used R data, as specified in the table. 

In the case of Johnson R-band data we assume F=R as 
discussed in Lugger (1989). For data in the Thuan and Gunn 
(1976) r and Johnson V we assume the average color of bright 
cluster galaxies implying F = r — 0.58 and F = V — 0.76 , 
as given by Schneider et. al. (1983) from which we also take 
the relevant k-corrections. Magnitudes have been corrected for 
the interstellar extinction Af by adopting the relation Af = 
0.07(csc 6—1) (Oemler 1974). Overall, the estimated magnitude 
error due to the uncertainty of the zero point of each plate plus 
the internal error is less than 0.2 mag . 



3. Luminosity function determination 

To compare the luminosity functions of different clusters it is 
important to define the galaxy samples of individual clusters 
according to uniform criteria, as discussed in L86. In partic- 
ular, the distribution of galactic types and luminosities varies 
from the cluster core to the field (Oemler 1974, Lugger 1986). 
Thus the luminosity functions were determined inside circular 
regions with a fixed radius of Ra = 1.7/z arcmin correspond- 
ing to 3 Mpc for H o =50 Km s~ Mpc~ , q = 1. Outside these 
regions a local field density was computed and compared with 
the field galaxy counts of Butcher fc Oemler (1985). The agree- 
ment was within 10 percent in 23 cases. Though the local de- 
termination is probably more appropriate, for the statistical 
comparison of different clusters we preferred to adopt the av- 
erage background to treat all the clusters in a uniform way 
and to obtain a closer comparison with the work of previous 
authors. In particular the background counts as a function of 
the red apparent magnitude mn were deduced from Butcher 
fc Oemler (1985) assuming R=F and an average color index 
J-F=1.0 as in Lugger (1989). A a straight line fit in the range 
14 < m R < 18 gives log N b = 0.503 • m R - 7.49, where N b is 
the number of background galaxies per square degree and per 
0.25 magnitudes interval. 

The galaxy samples of each cluster were corrected statis- 
tically by eliminating, in each magnitude bin of 0.2 mag, a 
number of galaxies estimated from the field density in the same 
magnitude bin. Finally a magnitude limit as close as possible to 
7713 + 3 was adopted in most (24) cases while for the remaining 
12 clusters, whose fainter magnitudes bins have an anomalously 
low population, a brighter limit was adopted (these objects are 
marked with # in column 8 of Table 1). 

The LFs were then fitted with a Schechter (1976) function: 

H L )d L = **(j F )\* P (-j F ) d (j F ) (1) 

maximizing the likelihood, as in Schechter &; Press (1976), by 
means of the MINUIT package of the CERN library. Prelim- 
inary results concerning a non parametric comparison of the 
luminosity functions have been presented elsewhere (Trevese 
et al. 1996). 

Since our aim is to analyze the relation between M* and 
the magnitude Mi of the brighter cluster member, each LF 
has been fitted with M* as free parameter and a fixed to the 
universal value —1.25 (Schechter 1976) and the first ranked 
galaxy has been excluded from the fit (see OH89). 
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Fig. 1. The luminosity functions in the F band of the 36 clusters of the sample. Each point represents the logarithm of the number n of 
galaxies in bins of 0.2 mag. The continuous curves represent the maximum likelihood fit with a Schechter function of the unbinned data. 
The brightest galaxy is not included in the fit and is not shown. 
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The resulting M* values for the entire sample of 36 clusters 
are reported in column 5 of Table 1 together with the Abell 
catalogue number, redshift (Struble fc Rood 1991) , Mi, Mio, 
Rood fc Sastry and Bautz fc Morgan classes , and the number 
of galaxies of each cluster used in the fitting. The source of 
morphological classification are Struble fc Rood (1987) for the 
RS type and Abell, Corwin fc Olowin (1989) for BM types. The 
differential LFs are shown in Figure 1 with the fitting function. 
It is to be noted however that the curve shown is not a fit to 
the points in the figure, which derive from an arbitrary binning 
of the data, but represents the maximum likelihood fit to the 
unbinned data. 

4. The statistics of Mi and M* 

The values of M* show a nearly gaussian distribution with M* 
= -22.66 ± 0.52 in agreement with the values -22.52 ± 0.45, 
-22.64 ± 0.50, -22.85 ± 0.23, found respectively by Dressier 
(1978), L86 and OH89. 



of the sample (see the discussion of L86 about the role played 
by A665 which has a somewhat uncertain value of M*). 

A straightforward interpretation would be that the selec- 
tive depletion of the bright end of the luminosity function, 
predicted by the galactic cannibalism model of Hausman and 
Ostriker (1978), does not agree with the observations. 

However, before deriving any physical conclusion from the 
observed correlation, it must be considered that the error on 
the zero point of the magnitude scale affects by the same 
amount both Mi and M* of the same cluster, causing a pos- 
itive Mi - M* correlation, even in the case the true values of 
Mi and M* are intrinsically uncorrelated. It is easy to show 
that, in the absence of any intrinsic Mi-M* correlation, the 
observed correlation coefficient would be r = <t c /(<ti • <t*), 
where <ti and a* are the observed standard deviations of Mi 
and M* respectively, and <r c is the standard deviation of the 
calibration error (see Massaro and Trevese (1996)). Thus, in 
our case a r.m.s. uncertainty e.g. as large as 0.3 mag on the 
calibration could account for the observed positive correlation. 
However, in a few cases it has been possible to find in the 
literature other photometric data on the galaxies we used to 
establish the zero point, obtained by different authors, some- 
times in different bands. A comparison of these data shows 
that the uncertainty is less than 0.2 mag. Thus it is possible 
that the calibration uncertainty accounts for the positive Mi- 
M* correlation only in part. In any case a reliable correction of 
this statistical bias would require an accurate estimate of both 
the random photometric noise and the calibration errors. 

A possible approach consists in studying magnitude differ- 
ences like (M* — Mio) and (Mi — Mio), which are independent 
of calibration errors, and of any global shift of the luminosity 
function. The choice of the 10-th ranked galaxy is motivated 
by the fact that the r.m.s. deviation of Mk has a minimum for 
k = 10 in our sample of 36 clusters. 



Fig. 2. M* versus Mi for the 36 clusters of the sample. The differ- 
ent symbols refer to the RS classes, as indicated. A positive global 
correlation is seen. 

Figure 2 shows the distribution of the 36 clusters of our 
sample in the Mi-M* plane. A fit of a straight line M* = 
a + aiMi to the data, assuming equal errors on both axes, 
gives oi = 1.12, or oi = 0.58 considering only the errors on 
M* . The correlation is positive with a coefficient r = 0.54 
and an associated probability P(> r) = 7 • 10~ of the null 
hypothesis. The clusters appear segregated according to their 
RS class, with early types towards the top left and late types 
towards the bottom right. 

Notice that L86, from a sample of 9 clusters, obtains oi = 
0.67 ± 0.37, r = 0.56 corresponding to P(> r) = 12%, while 
from the data of OH89, relative to 12 clusters, it is possible to 
derive oi = 0.22 ± 0.14, r = 0.45 corresponding to P(> r) = 
14%. Thus our result gives a statistically significant proof of 
the trend suggested by the previous findings of L86 and OH89. 
These results contrast with the findings of Dressier (1978) who 
obtains, for 12 clusters, oi ~ —0.5, r ~ —0.5, but a much 
higher probability of the null hypothesis due to the smaller size 



Fig. 3. (M* - Mio) versus (Mi - M w ) for the 36 clusters of the 
sample. The different symbols refer to the RS classes, as indicated. 
A negative correlation is seen. 

(M* - Mio) is plotted versus (Mi - Mio) in Figure 3 and 
appears negatively correlated: r = —0.46, P(~> \r\) = 4.6-10 -3 . 
The clusters are still segregated according to their RS classes. 
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Alternatively it is possible to perform a partial correlation 
analysis in order to single out the intrinsic correlation between 
Mi and M* , removing the correlation produced by a global 
shift of the luminosity scale. 

In the case of three (or more) stochastic variables Xi , i = 
1,2,3 , from the ordinary (zero order) correlation coefficients 
Tij ,it is possible to compute the partial correlation coefficients 
T iJ\h, (i,j,k = 1,2,3; i / j / k) defined by: 



)•(! 



(2) 



which represents the intrinsic correlation between Xi and xj 
corrected for the effect induced by the correlation of both Xi 
and Xj with Xk (see e.g. Anderson 1984). 

The ordinary correlation coefficients between Mi , M* and 
Mio are all positive : ri } *=0.54 , 7 , i,io=0.76 and 7 , »,io=0.89. 
The resulting partial correlation coefficient is ri }t -io = —0.50 
with an associated probability P(> \r\) = 1.9 • 10 -3 . 

The effect is statistically significant, thus providing a new 
constraint for any model of cluster formation and evolution. 

As already pointed out, in Figures 2 and 3 the clusters are 
segregated according to their RS class. The same effect appears 
using the BM classes. Thus we have divided the clusters into 
groups, corresponding to the RS classes F+I, C+L, B and cD 
respectively, to collect enough objects in each group. Then, for 
comparison, we have also grouped the clusters according to 
their BM class. 



Fig. 4. The average values < Mi > (lower panels) and < M* > 
(upper panels), computed for the subsamples corresponding to the 
Rood &; Sastry classes cD, B, C+L and F+I (left) and to the Bautz 
&; Morgan types (right). From late to early types < Mi > becomes 
brighter, while < M* > becomes cluster fainter. 

The average values < Mi > and < M* >, taken over each 
group, result negatively correlated. The effect is better seen in 
Figure 4 , where < Mi > and < M* > of the different groups 
are shown. It appears that < Mi > becomes brighter, while 
< M* > becomes fainter, in going from irregular to compact 
and regular clusters. 

A more precise characterization of this trend is obtained 
by computing, for the 36 clusters of our sample, the Kendall 



rank correlation coefficients r of Mi, M* and Mio with the RS 
and BM type, together with the relevant probabilities P(t) of 
the null hypothesis. The results are summarized in Table 2. 

A negative correlation appears of Mio with RS and BM 
types, but it is very weak and not statistically significant. It is 
possible to check from data in Table 1, that an error in the color 
transformation of 0.2 mag , e.g. causing a systematic shift of 
1/4 of the sample, does not appreciably change the correlation 
coefficients and the relevant probabilities. 

We also repeated the statistical calculations excluding from 
the sample the 12 clusters whose LFs have been determined 
in smaller range of magnitude. The results are essentially un- 
changed (in some cases the significance of correlation in even 
slightly improved). 

As a result we can conclude that we obtained for the first 
time significant evidence of a negative correlation of M* with 
both the RS and BM classes. We find also that the correlations 
oi M* -Mio and Mi-Mio with the morphological classes are the 
most significant. Finally we find a trend of Mi with the RS 
cluster type and we confirm the similar trend of Mi with the 
BM cluster type which was already known (see e.g. Sandage fc 
Hardy 1973). 

5. Conclusions 

We can summarize our results as follows: 

- The mean of the characteristic magnitudes M* , determined 

by maximum likelihood fits, is in good agreement with the 
values found in the literature. 

- M* is positively correlated with Mi with high statistical sig- 

nificance. 

- Clusters appear segregated in the Mi-M* plane according to 

both their Rood fc Sastry and Bautz fc Morgan type. 

- Including in the study also the magnitude Mio of the 10th 

brightest member, a partial correlation analysis shows a 
negative intrinsic correlation between Mi and M* . 

- We find statistically significant evidence that on average the 

magnitude Mi of the brightest cluster member is brighter 
in clusters of the earlier Rood fc Sastry and Bautz fc Mor- 
gan types. 

- The characteristic magnitude M* is on average fainter in 

clusters of the earlier Rood fc Sastry and Bautz fc Morgan 
types. The effect is statistically significant, providing a new 
constraint for theories of cluster formation and evolution. 

Once the RS types are interpreted as an evolutionary se- 
quence going from late type, irregular clusters to the more con- 
centrated, dynamically evolved cD clusters, the above results 
may support the cannibalism model of Hausman and Ostriker 
(1978). In this scheme, during the cluster evolution the first 
ranked galaxies grow by merging, becoming brighter. 

We stress that merging itself is not sufficient to explain 
the increase of M* . Rather one must assume that the merg- 
ing affects preferentially the most massive and bright cluster 
members. This also causes the increase of Mio, which however 
is found to be smaller. It is also important to remember that 
Mi is excluded from the maximum likelihood fit, as done by 
previous investigators, so that the determination of M* is not 
affected by the brightening of Mi along the RS sequence, from 
I to cD types, but simply measures the depletion of the bright 
end of the luminosity function. 

The emerging scenario could be the following. The positive 
global Mi-M* correlation which appears in the data is mostly 
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Table 2. Rank correlation statistics 
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intrinsic. Possibly because clusters are born with different lu- 
minosity functions which, to a first approximation, are glob- 
ally shifted towards brighter or fainter luminosities, according 
to some statistical distribution. During the subsequent evolu- 
tion the brightest member grows at the expense of other bright 
galaxies, which are most affected by dynamical friction, caus- 
ing the depletion of the bright end of the LF. This happens 
for both intrinsically brighter and fainter clusters. Once the 
spread in luminosity is reduced by computing the average over 
individual classes, the negative Mi-M* correlation appears. 

We stress that, beyond any evolutionary interpretation, our 
results show that the cluster LFs are not universal but depend 
on the cluster type in a systematic way. 

Intrinsic differences and evolutionary changes of the shape 
of the LF will be better understood through a non parametric 
analysis. 
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